Dissipation and Decoherence in a Quantum Register 
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A model for a quantum register TZ made of replicas of a d-dimensional quantum system 
(cell) coupled with the environment, is studied by means of a Born-Markov Master Equation (ME). 
Dissipation and decoherence are discussed in various cases in which a sub-decoherent enconding can 
be rigorously found. For the qubit case (d = 2) we have solved, for small A'^, the ME by numerical 
direct integration and studied, as a function of the coherence length of the bath, fidelity and 
decoherence rates of states of the register. For large enough the singlet states of the global 
au{2) pseudo-spin algebra of the register (noiseless at = oo) are shown to have a much smaller 
deoherence rates than the rest of the Hilbert space. 
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I. INTRODUCTION 



Preserving coherence in a quantum system is one of 
the most demanding features required to be able to take 
practical advantage from the implementation of the ob- 
jects of quantum information e quantum computation 
theory Indeed all the additional power, with respect 
to the classical case, arising from the quantum nature of 
the information-processing device depends on the com- 
plex linear structure of the state space of a quantum sys- 
tem and on the invariance of such structure under (uni- 
tary) time-evolution. The system is therefore endowed 
with a massive intrinsic parallelism and the capability of 
exhibiting interference. Unfortunately all of this holds 
only for closed quantum systems. Real systems are un- 
avoidably coupled with the environment in which they 
are enbedded, hence they have to be considered as open 
systems, no matter how weak is the interaction. The 
relevant state manifold has now a convex structure 
the dynamics in general is no longer unitary and the in- 
terference patterns may disappear. This set of effects 
is known as the decoherence problem The protec- 
tion of quantum-encoded information against environ- 
mental noise has been, up-to-now, mainly addressed in 
the framework of the so-called error correcting codes . 
These are essentialy schemes to encode redundantly in- 
formation in such a way that it can be recovered also 
when (a few) 'errors' due to external sources have oc- 
cured. Such schemes are often based on suitable mea- 
surement protocols that have to be performed frequently 
enough to keep the error level within the scope of the 
given encoding. Of course this implies that quantum 
information-processing systems have to be coupled with 
a classical measurement apparatus: even leaving aside 



the obvious pratical difficulties, such necessity naturally 
leads, at least, to a severe slow-down of the computa- 
tional speed. More recently ||^( see also 0|) has been 
put forward the idea that, conceptually, a more efficient 
quantum state protection can be realized by enconding 
the information in subspaces that the (non-unitary) dy- 
namics makes intrinsecally more robust against the per- 
turbation due to the environment. Here the attitude is, 
in some sense, opposite to that at the basis of error cor- 
recting codes: now one aims to encoding states that can- 
not be easily corrupted rather than to look for states 
that can be easily corrected. In this approach one has 
to assume explicit models of system-environment inter- 
action and try to design the various ingredients in such 
a way that that the algebraic-dynamical structure of the 
global system gives rise to the stable subspaces one is 
looking for. Since the typical environment consists of 
infinitely many degrees of freedom a direct Hamiltonian 
approach to the problem, is not the most suitable ex- 
cept for some simplified situation [||. In this paper we 
address the problem of dynamically stable quantum en- 
conding within a Master Equation formalism that allows 
us to deal directly with the marginal dynamics of the 
computational degrees of freedom. The relevant infor- 
mation about the environment is contained in a few pa- 
rameters appearing in the Master equation itself. The 
system considered is the model of a quantum register: 
N replicas of a given finite-dimensional quantum system 
(the cell). If the the cell is two-dimensional one obtains a 
N-qubit register. The key feature for the existence of the 
sub-decoherent codes is the possibility of partioning the 
register in clusters (possibly coinciding with a single cell 
or with the whole register) such that the cells within each 
cluster are collectively perturbed by the environnent. It 
is the dynamical symmetry of the cluster that allows to 
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single out collective (entangled) states that, at least on a 
short time scale, are unaffected by the noise and there- 
fore evolve unitarily. This mechanism has a well-known 
counterpart in quantum optics given by the phenomenon 
of subradiance ^ 

The paper is organized as follows: insect. II we introduce 
the model, in sect. Ill are discussed the general features 
of both Master Equation and the sub-decoherent codes. 
The cases of purely dephasing and dissipative coupling 
with the environment are analysed respectively in sect. 
IV and V. Sect. VI contains some conclusions. 



II. THE MODEL 

We call a system TZ a quantum register with N d- 
cells, if TZ is composed by N replicas of a d-level sys- 
tem. The Hilbert space is given by Tin = ^f^iTi-i, where 
Hi = (i = 1, . . . , TV) is the single d-cell Hilbert space. 
In particular, if d = 2 one has a N-qubit register. The 
set of the states (density matrices) of TZ is 



STz = {pe End(7^7^) ■■ P>0, p 



tr^p 



!}• (1) 



S-ji is not a linear subspacc of End(?i7^) but a convex 
submanifold. The register is coupled with an uncontrol- 
lable environment B (from now on the bath). The time- 
evolution of the states of the closed system TZ + B,is gen- 
erated by a Hamiltonian of the form H = H-ji + Hb + Hx- 
We now discuss the structure of each of these terms. The 
bath will be considered as a single bosonic field, namely 
iJg — ujk b\ bk describes a set of non-interacting lin- 
ear oscillator (field modes). The self-hamiltonian Hui 
of TZ, is assumed, for the time being, to be the sum of 
single-cell Hamiltonians Hf (i.e. the register is an array 
of non- interacting cells). The register-bath interaction is 
given by sum of the bath-cell interactions 



E 
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{gk^cblA^+h. c). 



(2) 



Here the Af^s are single-cell operators whose action is 
non-trivial only on the i-th tensor factor of Ti-jz, repre- 
senting the various interaction channels through which 
the i-th cell can be coupled with the bath. Although this 
kind of situation can be suitably handled by resorting to 
the notion of dynamical algebra |10|, in the following we 
shall assume that there is just one dominant interaction, 
the corresponding operator being Ai. As it is well-known 
the generic effect of Hj on the marginal dynamics of TZ 
is to induce dissipation and decoherence. The first ef- 
fect, of course, consists in the irrevesible loss of register 
energy into the bath. Decoherence is a pure quantum 
effect consisting in the destruction of phase coherence 
of the register states: due to the entanglement with the 
bath the initial pure preparations of the register become 
in a very short time-scale mixed. The interplay between 
these two phenomena is strictly related to the the nature 



of the operators {Ai}. Now we make another simplifying 
assumption, supposing that the A^'s are eigenvectors of 
the adjoint action of H-jz, 



[Htz, A,]^~eA, 



(e e R+). 



(3) 



This means that if e > 0, the necessarily non-hermitian 
and traceless, AiS (A|'s ) are the destruction (creation) 
operators of elementary cell excitations of H-ji. Notice 
that the energy e does not depend on the cell index i, 
in that we are considering replicas of the same physi- 
cal system. If one considers the zero temperature case, 
in which only the bath-vacuum is involved, the effect of 
Hx, therefore will be that of letting the register relax to 
the the A-vacuum \Aq), {Ai \Ao) = 0,Vi) by exciting the 
bath modes. On the other hand, if e = the possibly 
hermitian AiS belong to a symmetry algebra of H-n, and 
no energy-exchange at all occurs: the effect of register- 
bath interaction is pure decoherence. A quantity that will 
play an essential role in the following is the bath coher- 
ence length ^c, which, in an Hamiltonian approach, can 
be defined as the spatial scale over which the coupling 
constants gki = gk{i), have a non-negligible variation; 
when — oo, the guiS no longer depend on the qubit 
index i. This limit will be referred to as the replica sym- 
metric point, in that for = oo, the dynamics becomes 
invariant under the action of the symmetric group Sn of 
the cell permutations and only the collective operators 
A = Ai are effectively coupled with the bath. This 
situation corresponds to the well-known Dicke limit of 
quantum optics pl 

To exemplify this situation let us consider the S,c — oo, 
limit with Ai = . In this case, as far as the cou- 
pling with the environment is concerned, the relevant 
register operators are S"" = X^il^i'^f' ("^ — '^)- Let 
Hfi = eS^ + H]^, where is a qubit-qubit interaction 
term, and suppose that the latter is su(2)-invariant (i.e. 
[H]zi S"'] = 0, (a = ±,z)); one has then the commuta- 
tion relation [Hqz, S^] = zteS*^. There follows that for 
large time the register relaxes to the lowest S'^-eigenstate 
allowed by the total spin conservation. If H]^ = 0, this 
amounts to a ground-state relaxation, whereas if e = 
and H]^ ^ 0, there is no energy loss. This example will 
be discussed with greater detail in sect. V. 



III. MASTER EQUATION 

The quantum dynamics of the system TZ-\- B is highly 
non trivial, and exact results are difficult to obtain. Nev- 
ertheless one is mostly interested in the register marginal 
dynamics (i.e. forgetting about the bath degrees of free- 
dom) in order to study stability against external noise 
of the information-coding states of TZ. This issue can be 
conveniently addressed in the framework of the Liouville- 
von Neumann equation for open systems, the so-called 
Master equation (ME) . Following the standard Born- 
Markov scheme where one traces out the bath degrees of 
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freedom, which is assumed to be in the state pg, one ob- 
tains a closed equation for the marginal density matrix 
of TZ, of the form 

p = L(p)^(*adi/^ + L)(p), (4) 

where as usual ad H{p) = [p, H], denotes the adjoint ac- 
tion of H. The superoperator L is called the Liouvillian. 
The action of the non-Hamiltonian (dissipative) part is 

N-l p(-) 

L(P) = E {^1;^ A, PA]-^ (4 A,p + p A\ A, ) 

p(+) 

+ t'^^ApA, ~ (A, A] p + pA., A])}, (5) 

where the r^j^'^'s are temperature dependent coupling 
constants containing all relevant information about the 
bath. They are respectively associated to the process 
of de-excitation and excitation of the qubit system. At 

T = 0. one has F,-^^ = 0. The renormalized Hamiltonian 

' ''J 

H!j^ = H-ji -f 5Hii, where, by introducing the Lamb-shift 
parameters A^^^ 

5H^^ = ^(^17^4 A, + A, A]). (6) 

At zero temperature the excitation terms A^^^ are van- 
ishing. Notice that these terms make the cells interacting 
even though H-jz, is a free-cell Hamiltonian. On the other 
hand, it follows from relation (||), that [SH-ji^ H-ji] — 0; 
this means that, in this model, the Lamb shift terms are 
not responsible for additional register energy loss, but 
they are a source of dephasing. Let = tr''(6'l' 6^ pg) 
be the mean occupation number of the mode k in the 
initial (thermal) bath state ps- The explicit form for the 
coefficients appearing in the ME (P is 



k 




6 is the customary Heaviside function, and V denotes 
the principal part. From these relations it follows that 
r'^) and A'^^'^^ are hermitian. Furthermore F^^) > 
and r^") > r(+^. It is important to notice that the as- 
sumption (^) plays an essential role in the derivation of 
the ME, in that it allows to move to the interaction pic- 
ture (respect to Hti) Ai AiC"**^*. This is necessary 
in order to separate the (fast) dynamics generated by 
the self-Hamiltonian from the (slow) one generated by 
the coupling with the bath. When only the collective 
cell-operators A are coupled with the bath, H-ji has to 
satisfy condition (|^) only with respect to them. Given 
such H-ji one can obtain a family of new register Hamil- 
tonians, fulfilling the same constraint simply by adding 



terms commuting with {A, A^}. Introducing the nota- 
tion Af = 9{—a) Ai + 6{(t) A], equation (^ can be cast 
in the compact form 

Up) = ^ E (2 AfpA-^ - {AJ^ , p}) . (8) 

Diagonalizing the hermitian matrices F^""^ = 
llTij'!!, (c ~ ±) one obtain the following canonical form 
for the dissipative part of the Liouvillian |Q 

^(p) = l E A^([i>,ir] + [^^,pir])' (9) 

/j.cr— ± 

where {AJ^} are the eigenvalues of r*^°'). Moreover L'^ = 
^■u^ A^, denoting the components of the eigenvec- 
tors of r'^'^). The LZ^s will be referred to as the Lindblad 
operators. Given an initial pure preparation \ipo) of the 
register, one defines F{t) = (i/'o|p(i)|V'o) fidelity. Such 
quantity measures the degree of similarity with the ini- 
tial preparation that a state maintains during its time- 
evolution. Another quantity that one introduces in order 
to study the quantum coherence loss due to the bath is 
S{t) = tr {p{t) — p{t)'^), called linear entropy (or idempo- 
tency deficit). This quantity shares with the von Neu- 
mann entropy S = —trp log p, the fundamental property 
d[p] = <^ p^ = p (i.e. they both vanish iff p is a pure 
state). On the other hand, since the linear entropy does 
not involve trascendent operatorial functions, it is much 
simpler to evaluate than S. To characterize the degree of 
stability of the states it is useful to consider the short- 
times expansion 

oo 

<5(t) = 5(0) + E -(-)"• (10) 

In the following t„ {t~^ ) will be referred to as the n- 
th order decoherence time (rate). One straightforwardly 
finds 

l/r„" = -tr^ |e (fc) L"-^(p)L^(p)| {n > 1). 

(11) 

Since in the following ti will be play a major role, here 
we report explicitly the first decoherence rate 

l = -2tr{pL(p)}. (12) 

n 

In particular, for a pure initial preparation p = 
one has (5(0) = and 

/j.,cr— ± 

(13) 
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whereby one notices that the Hamiltonian component of 
the LiouviUian does not contribute to the first order deco- 
herence time (this comes from ti^{pa,dH^{p)} = 0). Of 
course these expressions obtained within the ME equa- 
tion formalism (which rehes on the Born-Markov assump- 
tion) differ from the ones that one could get by the exact 
temporal evolution induced by the interaction Hamilto- 
nian (H) (see for example [Q). Nevertheless, as far 
as the issue of code stability classification is concerned, 
this is not crucial in that the (exact) first order decoher- 
ence rate 1 / rf ^ is vanishing for pure initial state vanish- 
ing and essentialy corresponds to ti. 



A. codes 

The ME with initial condition p has formal solution 
p{t) — e*^{p), obtained by exponentiation of the Liou- 
ville super-operator L. The stationary solutions p{t) — p, 
are therefore the states belonging to kerL, where kcrL = 
{p e End{Hn) ■ L(p) = 0}. When t{p) = it follows, 
from equation (p^, that S{t) ~ 0{t ) (whereas for the 
fidchty one finds F{t) = 1 — Oit^)). Such a state will 
be called sub-decoherent. In general the adjoint action 
of H-ji maps sub-decoherent states onto states such that 
L(p) 7^ 0; but when S-ji n kerL is ad iJ^-invariant the 
LiouviUian evolution of each state p £ ker L becomes 
unitary: p{t) — exp(— i t) p exp(i H'^ t). In particular 
one has 5{t) = 0, Vi > (i.e. r^i = 0, Vn). This kind 
of state will be called noiseless. A subspace C C H-r, 
such that each density matrix over it is a sub-decoherent 
(noiseless) state will be referred to as a sub-decoherent 
(noiseless) code. 

Let us suppose to be sub-decoherent. First at all we 
notice that due to non-negativity of matrices F^"'^ and 
from the Schwartz inequality, each term of the sum in 
equation (|l3|) is non-negative. Therefore from — 
it follows that \\L1^ \^}f = IV')^, {^p,a), which 

in turn implies to be a simultaneous eigenvector of 
all the Lindblad operators. Conversely if is a simul- 
taneous eigenvector of the ij^'s then the sub-decoherence 

constraint t^^ = is trivially fulfilled. In other words 
a necessary and sufficient condition for the existence of 
a sub-decoherent code is the existence of a simultaneous 
eigenspace of all Lindblad operators LJ^ 

C = {IV^} e : IV^) ^ IV^), yp, a}. (14) 

The greater is d[a\ = dimCa the more efficient is the 
encoding. It is obvious that one has two quite different 
situations, depending on whether or not the ij^'s are her- 
mitian. In fact, if LJ^ = (ij^)^ one has that the Lindblad 
operators commute, [i^, LI] = J^ij'^i'^'j ^j] = 0; 
then there exists a non trivial Cq., furthermore Ti-ji = 
®aCa. On the other hand, if LJ^ 7^ (Lp'f the Lindblad 
operators no longer span an abelian algebra and cannot 
to be simultaneously diagonalized. The only candidate as 



sub-decoherent code is C = n^crkeri^. Indeed the Lind- 
blad operators satisfy relation (^), from which one de- 
rives that the only allowed eigenvalue is AJ^ = 0. The 
proof is as follows let {\Ei)}^-^ a iJ-Tj-eigenstates 

basis of Hn [Hn \E^) = E, \E,), E,+i > E„ D = d^). 
bmce the L+ are raising operators over the spectrum of 
H-jz one has L+ \Ei) oc where i' > i, Eii = Ei + e 

(in particular the maximum eigenvalue vector lEo) is 
annihilated by L+, (V/x)). Let — J2i=i^i\Ei) an 
eigenvector of with eigenvalue A 7^ 0; then one must 

have L+\ip) = Y.f^i' ^^^i = ^J2f=i \Ei) hence 
ci = 0. Acting on 1-0) with increasing powers of one 
analougously finds C2 = C3 = . . . ^ cd = 0, therefore if 
A 7^ one would have = 0. 

Let C the Lie algebra generated by the ij^'s (i.e. the 
minimal subspace of operators closed under commuta- 
tion containing {L'^}^a-) then the code C is nothing but 
the singlet sector of £; each G C is a one-dimensional 
representation space of C. From the general form of the 
Lindablad operators one has £ C (SfLi^i where Ci is 
the (local) Lie algebra generated by the Af^s. Generi- 
cally one has Ci = si (d, C), therefore if the above in- 
clusion is not strict it follows that C = {0}, which has 
no use for quantum encoding. In order to obtain mean- 
ingful codes one has to impose constraints on the alge- 
braic structure generated by the Lindblad operators. The 
smaller C the easier will be the task of finding (by repre- 
sentation theory) non trivial C. Notice that, given such a 
sub-decoherent code, if H!j^ belongs to the universal en- 
veloping algebra U{C) then C is also necessarily noiseless. 
The matrices T^'^'' and A^'^'' encode all the information 
about the spatial correlations among the register cells 
induced by coupling with the bath. The actual form of 
these correlations depends (see equation (|^)) on the de- 
tailed form of the coupling functions gki, on the bath 
density of the state and on temperature as well. 
Leaving aside strongly model-dependent considerations 
and in view of keeping the form of the ME here consid- 
ered as general as possible, in the following the matrices 
r(±)^ /\(±)^ be considered rather as a priori data of 
the problem defining the basic dynamical equation (Q). 
In other words they are treated as parameters that have 
to be 'engineered' in order to realize an advantageous sit- 
uation for quantum-enconding. In this context the bath 
coherence length is better defined in relation to the 
spatial behaviour of the F^^^ One can consider the fol- 
lowing particular regimes, corresponding to different 'ef- 
fective' topologies of TZ. It is just from these topologies 
that constraints on the algebraic structure arise. 

i) Tij — rSij, (V«,j) : this is the cell limit; the 
decoherence process occurs independently in each 
cell. The Lindblad operators coincide with the Af^a 
[C — (SfLiCi). 

ii) Tij — F (Vi, j) : this is the replica symmetric point; 
the decoherence is collective. The matrices T^'^'> 
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have constant entries, the only non-zero eigenvalue 
is A'^ and the corresponding Lindblad operators are 
given by L± = iV-1/2 Af (C^C,). 

The limit i) is the one usually considered in Error Correc- 
tion literature. The case ii) corresponds to the so called 
Dicke limit of quantum optics 

An interesting intermediate case between i) and ii) is 
when the register is partioned in clusters such that in 
each cluster the cells are coupled in the same way with 
the environment and different clusters are far enough to 
feel scorrelated environments. In other terms, if I (L) is 
the typical intra-cluster (inter-cluster) distance, we are 
supposing I <^ S,c {L £_c)- More formally we assume 
that there exists a partition {C\}fLi of the cell index set 
Njv, such that 

iii) Tij = Fq if i,j e C\, otherwise. The Lindblad op- 

erators are the cluster-ones L'^ = Ny^ SieCA ^ 
being N\ the number cells in the A-th cluster 
{C ^ ®fUCi). When M = N and M = 1 we 
recover respectively the cases i) and ii). 

For a clustered register the dynamics is invariant under 
the action of the group Q = Smi x • • • x ^ ^n', at 

the replica symmetric point (cell limit) one has G — Sn 
(G — {1}). Some comments are now in order. When 
the relation (||) holds we see that both in the hermi- 
tian and in the non-hermitian case the self Hamiltonian 
leaves invariant the code; nevertheless also in this rather 
special situation, due to the renormalizing terms (^) sub- 
decoherence does not necessarily imply noiselessnes. The 
point is that the F^^^'s and the A^^^'s in general cannot 
be diagonalized simultaneously. This can be understood, 
for example, by looking at the explicit form (^: in the 

matrix elements A^^' appears a sum over all the bath 

modes whereas in the F -^'^'s are involved only the modes 
degenerate with the single cell eigenvalue e. On the other 
hand, we see, from equation (^, that the leading con- 
tribution to A,-5^' comes from the same bath modes in- 

volved in F;if\ therefore assuming that A^f ^ and T^f^ 
ij ' ^ *j I'j 

have the same structure can be in many cases a good 
approximation. When this is the case also SH-ji can be 
written in terms of the Lindblad operators, namely each 
sub-decoherent code Cq is necessarily noiseless. 



IV. DECOHERENT COUPLING 

In this section we consider, the case in which the sin- 
gle cell-operators Ai in equation ^ are hermitian. Al- 
though this case is essentialy well-known we think that 
it is worthwhile to analyse it in that its exact solvabil- 
ity allows us to shed some light onto the general fea- 
tures of the decoherence process of many replicas of a 
given system coupled with the same environment. Here 



the ME is considered the starting point of the analy- 
sis, we do not assume any a priori relation like equa- 
tion (^. For the time being we set Htz = 0. Let 
\a) = |q;i, . . . , ctAf) denote a simultaneous eigenvector of 
the Ai's with Ai \a) = ai \a), {i = 1, . . . , N). The opera- 
tors I a) (a' I are eigenvectors of the Liouvillian 

L(|a)(a'|) = W{a,a') \a){a'\, 
Wia, a') = i (llalli - \\a'\\l) - ||a - a'\\l, (15) 

where ||/3||2, = (/3, M/3), (M = A, F = E.=± T^"^ e 
End(C^),/3 e C^). Notice that || • ||m, is a semi- 
norm only if M > 0, and a norm only if kerF = 
{0}. Each state over Hn can be written in the form 
p = J2a a' ^a,a' therefore the general solution 

of equation (||) is 

p(t) = ^i?„,a-e^("^"')*|a)(a'|, (16) 

a, a' 

whereby one derives the following expressions for fidelity 
and linear entropy 

F(i) = ^|i?„.fe^("-"')*, 

aa' 

<5(i) = l-;^|i?„„f e2**^("'"')*. (17) 

aa' 

By equation ( p^ the set of sub-decoherent and noiseless 
solutions of Liouville equation (||) is obviously related to 
the properties of the matrices A and F. First at all notice 
that from the second of equations (0), the immaginary 
terms in equation ( p^ play no role in decoherence (in the 
restricted meaning): indeed they give rise to the unitary 
transformation 

C/A(i) = e-'**^^ =^e'll"ll^*|a)(a|. (18) 

a 

It is straightforward to verify that the linear entropy is a 
monotonic non-decreasing function of time, indeed 

m^2j2\Ro.a'\^\\a-a'\\l>0, (19) 

aa' 

the inequality following from the non-negativity of F. 
Cases i) and ii) imply, from W{a, a) — 0, that the di- 
agonal states Pa = are fixed points of the Liou- 
villian evolution. Furthermore if a — a' € kerF, one 
has that the real part of W{a, a') vanishes. Case i) cor- 
responds to a solution that one could obtain assuming 
that each cell is interacting with its own independent 
environment. From equation ( |l5| ) it follows that the 
maximum decay rate is 0{N). In case i) kerF — {0} 
and only a — a' survives. If the single-cell eigenvalues 
ai are non-degenerate the eigenspace Hdai}) is one- 
dimensional and therefore useless for quantum encod- 
ing. If instead the ai's are TOi-fold degenerate, then 
d[a] = dim7i({ai}) = JliLi ™i- The density matrix cor- 
responding to IV') G H{{ai}) evolves according the uni- 
tary transformation UA{t)- these states are noiseless. 
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The largest dimension for the noiseless code W({ai}) is 
obtained for ai = Um, (Vi) where aM is the single-cell 
eigenvalue with the maximum degeneracy. In the qubit 
case Ai = af and ai = ±1/2. In case even A^^^ is pro- 
portional to the unit matrix, then the unitary transfor- 
mation (|l^ ) becomes trivial, being X]i Q^f = ^/4, (Va). 
For the initial state jV'o) = 2"^^^ |cr), uniform lin- 
ear superposition of all the basis states, one can obtain 
explicit analytical expressions for the linear entropy and 
the fidelity 

5(t) = l_e-r^*cosh^(ri), 
F{t) = cosh^(r/20. (20) 

For t ^ oo one finds F ~ 2"^ and ^ - 1 - 2"^, re- 
sults that can be immediately understood from p{oo) ~ 
2-N \a){a\. Let us turn to the case ii). The operator 
A = J2i plays the role of pointer observable ||l^: the 
diagonal elements with respect to its eigenstates basis 
of the density matrix do not decohere, whereas the off- 
diagonal decay with a rate that is proportional to their 
distances from the diagonal. Now dimkerF = — 1, 
and the no-damping condition becomes J^i'^i — J^i'^i- 
This means that in that case the space Tia spanned by 
the set Ba = {!«) : J2f=i '^i — decoherence-free. In 
passing we note that, since A is an extensive observable, 
at the replica-symmetric point the maximum decay rate 
is 0{N'^). 

In case iii) the matrix T is block constant and a — 
a' e kerr, iff E.^c, « = E.ec. = 1, ■ • • , M). 

Now the relevant operators are the cluster operators 
Lx — X^jeCA states built over a simultaneous 

eigenspace of the L\s evolve in a noiseless way. As usual, 
the situation is best exemplified by the qubit case. Let 
us assume that Ai = erf, and N even. At the = oo 
point the most efficient noiseless-enconding is obtained 
building states over the eigenspace — Q. If F is par- 
tioned in blocks of m (even) elements one can encode 
in the subspace with zero cluster z-spin. Such code has 
dimension 

dm . (21) 

This enconding, with m = 2 is essentialy that proposed in 
I p^ . Till now we have supposed that the self-Hamiltonian 
were vanishing. If this is not the case, one has that for 
an initial noiseless preparation the state evolves infinites- 
imally in a unitary fashion. For finite time the (possible) 
non-commutativity between H-ji and the relevant Lind- 
blad (cell, cluster, register) operators, destroys the coher- 
ence of p. When relation (^) holds (e = 0) H-n commutes 
with the Lindbald operators. Working in a basis that si- 
multaneously diagonalizes Hn and the A^ 's one sees that 
U/^{t) — > exp(— z t Hti), the evolution will remain unitary 
for finite times; the initial pure states never get mixed. 



V. DISSIPATIVE COUPLING 

In this section we consider the case of non-hermitian 
Ai] namely the case when the relation (^) holds with 
e > 0. At zero temperature the eigenvalues are vanish- 
ing. On the other hand since F^"^ > 0, A~ > 0, (V^) one 
can immediately check that the register energy E-jiit) = 
tr^ {p{t) H-jz) is a monotonic non-increasing function. In- 
deed 

En{t)=tv'^ (t{p)Hn) 

= -e^A;tr^(i+L;p) <0, (22) 

where we have used the irrelevance of Hamiltonian com- 
ponent of L (that is ti^ {H-ji, [H!j^^ p]) = 0,) the rela- 
tion [ij^i^'^j H-ji] = 0, (that follows form equation (||) 
which holds for the Lindblad operators as well), and the 
non-negativity of operators and p. As it has been 

observed in sect. Ill in the present case a sub-decoherent 
code can be obtained if C = fl^keri^ ^ {0}. 
Restated in this formalism the essence of the result of 
reference |^ for the qubit case is that at the — oo 
point the Lindblad operators (and the renormalized self- 
Hamiltonian) belong to a iV-fold tensor representation 
of a semisimple (dynamical) Lie algebra, out of which a 
non-trivial C can be built when N is large enough. In 
the cell limit i) if one can find a subspace Ci C Ti. anni- 
hilated by both A^'^^ and A^'^ then C = Ci (g) • • • ® Cat. 
An analog construction can be made in the cluster limit. 
An important example is given by the qubit case. One 
can design a register that supports noiseless encondings 
if one is able to build TZ in such a way that iii) is satisfied 
with m = 4 qubits for cluster. Then, according reference 
to 1^, a logical qubit can be encoded in each cluster. It 
is important to note that the dimension of C decreases 
passing from ii) to iii), and from iii) to i). 

In general one has F^^^ = T^^\i, j). The first order 
time-scale ti is a functional of Itp), depending on ^. The 
optimal states, with respect to the storage reliability 
on short times, are those that minimize this functional 
for a given bath coherence length. Let us assume that 
Tij = ro74(* ~ i)' where 75 (cc) — > 1, when ^ 00 and 
7j(a;) — > Sx,o, when — >■ 0"*". The latter situation cor- 
responds to the case in which each cell is coupled with 
an independent bath, therefore G (0, 00) interpolates 
between the independent bath limit i) and the infinite 
coherence length bath case ii). 

A. qubit case 

Now we specialize to the d = 2 case: Af — af. Let 
the self-Hamiltonian be of the form Hn = eS"^ -I- 
where the second term is a qubit-qubit interaction. In 
quantum computation applications such a kind of term 
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might arise, for example, during the gate processing. If 
we assume that [H^, S"] = 0, (a = z,±) then equation 
holds. Now we briefly recall the result of reference H] 
at the replica symmetric point. When — oo, one finds 
that: 



i) the total spin operator S"^ = (5^)^ 
is a constant of the motion. 



l/2{5+, S-}, 



ii) defining in the obvious way a iSjv-action T over the 
density matrices manifold S-ji, one has Tg-LTj = 
L, (Va G Sn), 

iii) the Lie algebra C generated by the Lindblad oper- 
ators is nothing but the global su{2). 

iv) Since at = oo the coupling functions gki are as- 
sumed to be strictly qubit independent also the 
Lamb-shift matrices A,-if^ have constant entries 

Aq^ViiJ). The renormalizing term 



can then to be written as SH-jz — Ag 5*+ 5 



A+5- 



S+ 



From i)-iv) it follows that the Hilbert space TC-r, splits dy- 
namically according the Clebsch-Gordan decomposition 
of the n-fold tensor representation of su{2) 



^N{S) 

~=1 



UriS), 



(23) 



where Smin = ( Smin — 1/2) if N is even (odd). 
The subspace Hr{S) is an irreducible SM(2)-module cor- 
responding to the total spin eigenvalue S {S + 1), the 
latter occurring with multiplicity 



nN{S) 



i2s + i)m 



{N/2 + s + iy.{N/2-sy: 



(24) 



The general state over HriS) has the form: p = 
Em m'=-s Pm,m' \SM){SM'l where \SM) = S{S + 
1) |5M), \SM) = M \SM), [M = -S,...,S) and 
analogously for \SAI'). For a pure state one has 

n(oo)-i = 2r(-^ ((^1 s+s- \^) - m s- m') (25) 
+ 24+\{i^\S- S+\^)-MS+m'). (26) 

In particular, if = \SM) one obtains {2ti)^^ = 
ToCliS, M) + T^Q^^ CliS, M), where Cl{S, M) = 
5(5 + 1) - M(M ± 1). Let us consider the zero- 
temperature case (Fq^"* = 0) when only the de-excitation 

processes with strength proportional to Fq ^ are active. 
If is a lowest-weight spin state (i.e. 5^ \ip) = 0), 
one has ti(oo) = oo. This result is true for all the 
decoherence-times t„. At finite temperature the (exci- 
tations) terms weighted by Fq^'' are present as well. 

On the su{2)-singlets IV') € C ee ©"f f ^7i;,.(0), one has 
5+ \ip) — S~\tp) = 0, and SH-ji = 0; furthermore from 
the su(2)-invariance oi H-j^ there follows that the unitary 



part of L maps the singlet sector onto itself, namely C 
is noiseless. From equation (^^ it follows that the min- 
imum cluster size to encode a noiseless logical qubit is 
= 4. Defining (in obvious binary notation) the states 
\A) EE 10011) + |1100), \B) EE 10110) + |1001), |C) = 
|1010) + |0101), an orthonormal basis of C is given by 



|0) 

ID 



2- H\B)-\A)), 

3- 1/2 2-1 |B)). 



(27) 



If — this two states are energy degenerate, for non- 
vanishing qubit-qubit interaction the degeneracy is lifted. 
For example if 



J E i^^^s + <^J + <7+)}, (28) 



<ij> 



is a Heisenberg coupling between nearest neighbour 
qubits arranged on a ring topology, one finds that |0) 
and |1) are energy eigenstatcs with eigenvalues respec- 
tively given hy Eq = J and Ex — — J. Since is su{2)- 
invariant it is always possible to choose the singlet 
among its eigenvectors. It should be emphasized that the 
su{2) singlet sector is noiseless for a wider class of ME's, 
with Lindblad operators (and self-Hamiltonian) given by 
arbitrary functions of the global operators 5", (a = ±, z) 
p7| . Indeed if these operators have the form 



X = cxl + F{{S'^}) 



(29) 



where F is an arbitrary operator-valued analytic func- 
tion, then - since F I?/;) — 0, (V 1-0) G Cjv)- one obtains 
A|cjv = ci I. This latter condition is sufficient to preserve 
the sub-decoherence of Cat. Another way to understand 
this result is that the operators described by Eq. ( p9| ) co- 
incide with the iSAT-invariant sector (symmetic subspace) 
of End(7iK) ■ Since Cn is an irreducible iSjv-module from 
the Schur lemma it follows that A|cjv oc I. 
Turning back to the general case ^ G (0, oo), if A^ = 2, for 
the initial states | tt), I ii), \^t,s) - 2-1/^ (| ti) ± | it)), 
one immediately finds 



r/^(e) = (2F(-V\ r/^(e)-(2FW) 
r^-^(0 = (2(F[,-)+FW)(l±7c(l)))-i 



(30) 



These equations show that in the generic case (F^^) ^ 0) 
for finite coherence length ^ all the first order decoher- 
ence times are finite as well, whereas for ^ — > oo the 
singlet (^), diverges with ^. Of course for this latter 
state, since L^=oo(|V''*)('0'*l) — 0, all the t„'s diverge. 
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FIG. 1. First order time-scale for ttie symmetric 
state It/)"^™} = (S'+)2|0), and a singlet state, iV = 4 
r(j,j)=0.1e' 



In the general case when the matrices T'^^'' are not 
block constant one has to resort to numerical calcu- 
lations. We have solved equation by direct nu- 
merical integration in the qubit case with H-ji = eS^. 
Rather than using the form (0) for the ME parame- 
ters, we have choosen a phenomenological parametriza- 
tion such as r|^^ = rQ"*^^ ^-b-Jl/ic ^j^j neglected the self- 
Hamiltonian renormalization. In figure is reported 
the behaviour of ti(^c) for a = 4 singlet and the 
highest-weight SM(2)-vector belonging the S = 2 multi- 
plet. We see that for a wide range of the decoherence 
time of the singlet state is much larger than that of the 
symmetric state. 




FIG. 3. Linear entropy as a function of the time for the 
= singlet and triplet state (iV = 2, To = 0.1 e'l'~^l/«) 
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FIG. 4. Fidelity difference, between a A'^ = 4 singlet 
and the state (S^)'^ |0), for different bath coherence lengths ^ 
(Fo = 0.1e-l'-^l''«.) 
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FIG. 2. Fidelity as a function of the time for the = 
singlet and triplet state (iV = 2, Fq = 0.1 e"''"-''/^) 
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FIG. 5. Linear entropy difference, between the state 
(S"*")^ (0), and a A'' = 4 singlet, for different bath coherence 
lengths ^ (Fo = 0.1 e" I* 



In figures (^) is compared the behaviour of the 
fidelity and linear entropy of the N = 2 singlet and 
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{S"" = 0) triplet states at finite ^c- Figures (§) and (||) 
show, as a function of time, the difference of fidelity and 
linear entropy, between one of the = 4 singlets and the 
symmetric state (S*"*")^ |0) € 7ii(2), for various bath co- 
herence length. These simple calculations strongly sug- 
gest that the noiseless enconding at infinite coherence 
length remains, for sufEciently large ^c, more robust than 
all other states 



B. Gauge transformation 

We end this section by showing that for a class of non- 
trivial qubit couplings, connected to the limits ii) and iii) 
via a local gauge transformation, it is possible to build 
subspaces annihilated by the dissipative component of 
the Liouvillian. This transformation is a generalization 
of the one considered in Here we give the proof 

for the replica-symmetric case, the cluster case being 
a straightforward generalization. Let us suppose that 
r^(^=^^ = r(±) e*('^(*)-'^(J», where r(±), e R, <I):Nn ^ R. 
This kind of situation is not completely fictitious: for 
9ki ^ e**"'""' when there is just one bath mode k degen- 
erate with the qubit energy e, from the first of equations 
(0) follows that r''^'' - (,^k{r,~r,)_ Introducing the op- 
erators = X^jLi ^"^^"'''^j The dissipative Liouvillian 
has the canonical form (^) with {X^} — {T'^'^^} and Lind- 
blad operators given by the i^'s. The operators {L^j^, 
spanning a Lie Algebra isomorphic to A generated 
by {L'^jcr, are obtained from the latter by means of the 
(local) U{1) gauge transformation 

: Endinn) ^ End(7^7^): X^U^X 

N 

= exp{ze~i } e <=iC/(l)„ (31) 

where we recall that Hf is the single cell Hamiltonian 
fulfilling relation (|3|) with the ^i's. The unitary operator 
U e End(HK) maps the singlet sector C of ^ onto the 
one oi A^. Therefore p £ C ^ L{T^p) = 0. The new code 
U^{C), is noiseless depending on the transformation prop- 
erties of H'^ under T^. If H-jz — T^{Hti), (local gauge 
invariance) it follows that Utj,{C) is noiseless under iff 
C is noiseless under L^^o (replica independent case). Let 
us consider, for example, the qubit case with N — 2 and 
(f){j) — (f) j {(f) £ Ti) and H!j^ = e S^. The noiseless state is 
now the singlet \tps) = 2^^/^{\ tJ,) — | |f)). It is mapped 
by onto U^H^s) = 2-i/2(e'*/2| ||) _ e^^4'/2 \ ||))^ 
particular for (j) = n, one has t/jrlV's} = li^t), that is the 
triplet state becomes the noiseless one. For e (0, tt) one 
has a smooth interpolation from \ips) to {ipt)- It should 
be emphasized that even if T^{Hti) = H-jz, generally one 
has that the many-qubit correction 5Hti is not invariant. 
Nevertheless the Hamiltonian part of L, does not affect 
the first order decoherence rate: Ucf,{C) is sub-decoherent. 



VI. CONCLUSIONS 



In this paper we have studied a model of quantum reg- 
ister TZ with N cells made of replicas of a d-dimensional 
quantum system. The register TZ is coupled with the en- 
vironment, modelled by a thermal bath of harmonic os- 
cillators, through single-cell operators Ai . The latter are 
step operators over the spectrum of the cell Hamiltonian. 
The reduced dynamics of TZ is studied by a Master Equa- 
tion (ME) obtained in the Born-Markov approximation. 
The ME provides a very natural and powerful tool to dis- 
cuss, in a unified way, the various aspects of decoherence 
and dissipation phenomena induced in TZ by the bath. 
The effect of the environment splits into two contribu- 
tions: a renormalization of the register self-Hamiltonian, 
that makes the cells effectively interacting, and an irre- 
versible component describing the decay processes. The 
latter can be cast in canonical Lindblad form by diago- 
nalizing the N xN matrices r^^^^ ; which contain all infor- 
mation about the effective spatial structure of TZ in the 
given environment state. Three situations which appear 
to be relevant for quantum encoding have been discussed: 
i) all the cells are coupled with the environment in the 
same way, ii) different cells feel different environments, 
iii) the register can be decomposed in uncorrelated clus- 
ters, such that the cells within each cluster satisfy i). In 
each of these cases one can show the existence of sub- 
spaces C such that an initial pure preparation 1-0) G C 
has vanishing linear entropy production rate. The states 
in C therefore - on a short time-scale - maintain quantum 
coherence: C can be thought of as a subdecoherent code. 
The latter is obtained as simultaneous eigenspace C of the 
Lindblad operators L^, given by linear functions of the 
A's associated with the register cells. Depending on the 
structure of the Lie algebra C generated by the L^^'s one 
has to face rather different situations. For A hermitian 
C is abelian, the Hilbert space splits in a direct sum of 
the simultaneous eigenspaces C : the ME is exactly solv- 
able. Analytical expressions for decoherences rate can be 
found in the qubit case. In the non-hermitean case C 
is non abelian, the Hilbert space splits according to the 
£-irreps, C is the common null space of the i^'s (singlet 
sector of C). The latter exists, according to reference 
if the size of the clusters satisfying i) is large enough. 
For the qubit case the minimun cluster size required to 
encode one logical qubit is = 4 : a register made of 
M clusters of four qubit each supports a 2^^-dimensional 
subdecoherent space. If C is left invariant by the renor- 
malized self-Hamiltonian of TZ the time-evolution of 
the sub-decoherent states is unitary: the code is noise- 
less. In this case the relevant algebra is C generated by 
the Lindblad operators plus H'-ji- Furthermore we have 
shown that there exist cases with non-trivial cell depen- 
dence that can be mapped onto ii)-iii) by a suitable lo- 
cal gauge transformation. The degree of stability of the 
resulting codes depends on the covariance properties of 
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the renormalized self-Hamiltonian. When the F^^^'s are 
not block-diagonal one has to resort to numerical cal- 
culations. We have integrated the dissipative ME of a 
qubit register. The results show that for a wide range of 
bath coherence lengths the singlet states (noiseless at 
— oo) are more robust, namely their entropy increases 
more slowly on the time scale of decoherence. 
The problems related to the pratical realizations of the 
registers satisfying the constraints for the suggested en- 
codings, the preparation as well as the gate manipula- 
tions of the codewords necessary in the quantum compu- 
tation applications, are of course open issues that deserve 
further investigations. 
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